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Abstract. Boolean network models of molecular regulatory networks 
have been used successfully in computational systems biology. The 
Boolean functions that appear in published models tend to have spe- 
cial properties, in particular the property of being nested canalizing, 
a concept inspired by the concept of canalization in evolutionary biol- 
ogy. It has been shown that networks comprised of nested canalizing 
functions have dynamic properties that make them suitable for mod- 
eling molecular regulatory networks, namely a small number of (large) 
attractors, as well as relatively short limit cycles. 

This paper contains a detailed analysis of this class of functions, based 
on a novel normal form as polynomial functions over the Boolean field. 
The concept of layer is introduced that stratifies variables into different 
classes depending on their level of dominance. Using this layer concept 
a closed form formula is derived for the number of nested canalizing 
functions with a given number of variables. Additional metrics consid- 
ered include Hamming weight, the activity number of any variable, and 
the average sensitivity of the function. It is also shown that the average 
sensitivity of any nested canalizing function is between and 2. This 
provides a rationale for why nested canalizing functions are stable, since 
a random Boolean function in n variables has average sensitivity ^ . The 
paper also contains experimental evidence that the layer number is an 
important factor in network stability. 



1. Introduction 

Canalizing Boolean functions were introduced by S. Kauffman and collab- 
orators [20] as appropriate rules in Boolean network models of gene regula- 
tory networks. More recently, a subclass of these functions, so-called nested 
canalizing functions (NCF) was introduced [21 J and studied from the point 
of view of stability properties of network dynamics. A multi-state version 
of such functions has been introduced in |33[ I34j . where it was shown that 
networks whose dynamics are controlled by nested canalizing functions have 
similar stability properties, namely large attractor basins and short limit 
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cycles. An analysis of published models, both Boolean and multi-state, of 
molecular regulatory networks revealed that the large majority of regula- 
tory rules in them is canalizing, with most of these in fact nested canalizing 
[Ml [22l [Ml |33]. Thus, nested canalizing rules and the properties of net- 
works governed by them are important to study because of their relevance 
in systems biology. Furthermore, they are also important in computational 
science. In jl6J it was shown that the class of nested canalizing Boolean func- 
tions is identical to the class of so-called unate cascade Boolean functions, 
which has been studied extensively in engineering and computer science. It 
was shown, for instance, in [8] that this class has the property that it corre- 
sponds exactly to the class of Boolean functions with corresponding binary 
decision diagrams of shortest average path length. Thus, a more detailed 
mathematical study of nested canalizing functions might have applications 
to problems in engineering as well. 

In this paper, we carry out such a detailed study for the case of Boolean 
nested canalizing functions, obtaining a more explicit characterization than 
the one obtained in [TTj. We introduce a new concept, the layer number, 
leading to a finer classification and, in particular, an explicit formula for the 
number of nested canalizing functions. This provides a closed form solution 
to the recursive formula derived in p^Y] , which may be of independent math- 
ematical interest. We also study standard properties of Boolean functions, 
such as variable activity and average sensitivity. In particular, we obtain 
a formula for the average sensitivity of a nested canalizing function with n 
variables, and show that, for all n, it lies between 2!?=^ ^ lower bound and 
2 as upper bound, which is much smaller than ^, the average sensitivity 
of a random Boolean function in n variables. This can be interpreted as 
providing a theoretical justification for why Boolean networks with nested 
canalizing rules are stable. We also find a formula for the Hamming weight 
(the number of I's in its truth table) of a nested canalizing function. Finally, 
we conjecture that a nested canalizing function with n variables has maxi- 
mal average sensitivity if it has the maximal layer number n — 1. Based on 
this result, we conjecture the tight upper bound | for this value. The paper 
is organized as follows. We first review existing results on nested canaliz- 
ing functions and networks, after which we introduce some definitions and 
notation. The subsequent sections contain the main results of the paper. 

2. Background 

In j32j it was shown that the dynamics of a Boolean network which op- 
erates according to canalizing rules is robust with regard to small pertur- 
bations. In [19], an exact formula was derived for the number of Boolean 
canalizing functions. In [30], the definition of canalizing functions was gen- 
eralized to any finite field Fg, where (; is a power of a prime. Both exact 
formulas and asymptotic values for the number of generalized canalizing 
functions were obtained. 
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One important characteristic of (nested) canalizing functions is that they 
exhibit a stabihzing effect on the dynamics of a Boolean network. That 
is, small perturbations of an initial state do not grow larger over time and 
eventually end up in the same attractor as the initial state. This stability 
is typically measured using so-called Derrida plots which monitor the Ham- 
ming distance between a random initial state and its perturbed state as both 
evolve over time. If the Hamming distance decreases over time, the system 
is considered stable. The slope of the Derrida curve is used as a numerical 
measure of stability. Roughly speaking, the phase space of a stable system 
has few components and the limit cycle of each component is short. 

In [22], the authors studied the dynamics of nested canalizing Boolean 
networks over a variety of dependency graphs. That is, for a given random 
graph on n nodes, where the in-degree of each node is chosen at random 
between and k, for k < n, a nested canalizing function is assigned to 
each node in the in-degree variables of that node. The dynamics of these 
networks was then analyzed and the stability measured using Derrida plots. 
It was shown there that nested canalizing networks are remarkably stable 
regardless of the in-degree distribution and that the stability increases as 
the average number of inputs of each node increases. 

Most published molecular networks are given in the form of a wiring di- 
agram, or dependency graph, constructed from experiments and prior pub- 
lished knowledge. However, for most of the molecular species in the network, 
little knowledge, if any, could be deduced about their regulatory rules, for 
instance in the gene transcription networks in yeast [15] and E. Coli [3J. 
Each one of these networks contains more than 1000 genes. Kauffman et. al 
[2Tj investigated the effect of the topology of a sub-network of the yeast tran- 
scriptional network where many of the transcriptional rules are not known. 
They generated ensembles of different models where all models have the 
same dependency graph. Their heuristic results imply that the dynamics 
of those models which used only nested canalizing functions were far more 
stable than the randomly generated models. Since it is already established 
that the yeast transcriptional network is stable, this suggests that the un- 
known interaction rules are very likely nested canalizing functions. Recently, 
a transcriptional network of yeast, with 3459 genes as well as the transcrip- 
tional networks of E. Coli (1481 genes) and B. subtillis (840 genes) have 
been analyzed in a similar fashion, with similar findings [^. 

The notion of sensitivity was introduced in [12]. The sensitivity of a 
Boolean function of a variable x is defined as the number of Hamming neigh- 
bors of X on which the function value is different from that on x. The average 
sensitivity of the function is then computed by taking the average value of 
the sensitivities of the function on all possible input values x. Although the 
definition is straightforward, the sensitivity measure is understood only for 
a few classes of functions. In [39], asymptotic formulas for a random mono- 
tone Boolean function are derived. Recently, S. Zhang [l5] found a formula 
for the average sensitivity of any monotone Boolean function, and derived a 
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tight bound. In [lO], I. Shmulevich and S. A. Kauffman obtained the average 
sensitivity of Boolean functions with only one canalizing variable. In pT] , 
Layne et al. studied network stability of partially canalizing functions using 
the average sensitivity of such functions. 

3. Definitions and Notation 

To set the stage for this paper and for the sake of completeness we restate 
some well-known definitions; see, e.g., [21]. Let F = F2 be the field with 2 
elements, and let / : F" ^ F. It is well known [31] that / can be expressed 
as a polynomial, called the algebraic normal form (ANF) of /: 

/(xi, . . . , Xn) = CLki...kn-^l i . . . Xn 

0<fcj<l, 
i=l,...,n 

where each coefficient aki...k„ ^ ^- The number ki + k2 + • ■ ■ + kn is the 
multivariate degree of the term ak^k2...knXi^^ ^2^^ ■ • ■ Xn^" for each nonzero 
coefficient akj^k2...kn- The greatest degree of all the terms of / is called its 
algebraic degree, denoted by deg{f). The symbol © stands for addition mod- 
ulo 2, whereas the symbol -|- will be reserved for addition of real numbers. 
First we need a technical definition. 

Definition 3.1. The function f{xi,X2, ■ ■ ■ is essential in the variable 
Xi if there exist xj, . . . , x*_^ , x*^^, . . . , x* G F such that 

f{xi, . . . , x^_i, 0, x^_^_l, . . . , x„) 7^ fi^i^ • • • ) 1, x^_^_i, . . . , x„). 

The next two definitions state the requirements for Boolean functions to 
be canalizing, respectively nested canalizing. The concept of canalization in 
gene regulation goes back to work of the geneticist C. Waddington in the 
1940s [43], who developed it as a possible answer to the question of why 
the outcome of embryonal development leads to predictable phenotypes in 
the face of widely varying environmental conditions. Canalized traits of an 
organisms, those that are stable under (some) environmental perturbations, 
are phenotypically expressed only in certain environments or genetic back- 
grounds. The regulation of the genes responsible for the development of 
such traits by other genes has to be able to buffer these perturbations. In 
[23], S. Kauffman tried to capture the spirit of these features in the context 
of Boolean network models of gene regulatory networks. In that setting, a 
Boolean function / is canalizing in a variable x, with canalizing input a and 
canalized output b, if, whenever x takes on the value a, then / outputs the 
value b, regardless of the states of the other variables in /. Boolean network 
models built from functions with this property have been shown to have 
dynamic features that match those of gene regulatory networks. 

Definition 3.2. Let a,b G ¥. A function /(xi, X2, . . . , x„) is < i : a : b > 

canalizing i/ /(xi, . . . , Xi_i, a, Xj+i, . . . , x„) = b, for all xj, j ^ i, where 
ie {l,...,n}. 
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Thus, if the function receives its canahzing input a for variable Xj, then 
the function obtained by substituting a for Xi becomes a constant function 
equal to h. 

The motivation for the next definition is the general stability of gene 
regulatory networks. In the context of a Boolean representation of gene 
regulation, if a gene Xi does not receive its canalizing input a, then, in 
principle, the function obtained by substituting a© 1 for Xi can be a random 
Boolean function, with uncertain stability properties. In order to remedy 
this deficiency, it was proposed in [21] that in this case there should be 
another variable, Xj that is canalizing for a particular input; and so on. The 
next definition captures this intuition. 

Definition 3.3. Let f he a Boolean function in n variables. Let a he a 
permutation of the set {1, 2, . . . , n}. The function f is a nested canalizing 
function in the variable order a^o-(i) , . . . , x„(^n) with canalizing input values 
ai,...,a„ and canalized values bi,...,bn, if it can he represented in the 
form 



f [xi, . . . , Xn) — 







= ai, 




b2 


^-(1) 


= oT, 






^a(l) 


= al. 


Xa{2) = 02,3^(7(3) = ^3, 


bn 




= oT, 


• • • 1 '^(7(n— 1) ~ O^ri— li Xa{n) ~ 


[ bn 




= oT, 


■ ■ ■ 1 '^(7(n— 1) ~ Ojn—l-: Xa{n) ~ ^n- 



Here, a = a ® 1 . 

The function f is nested canalizing if f is nested canalizing in the variable 
order Xo-(i), . . . , ^^-{n) for some permutation a. 

Let a = (ai, a2, . . . , an) and /3 = (61, 62, • • • , bn)- We say that / is {a : q : 
/3} NCF if it is NCF in the variable order 

3^o-(i) ) • • • ) ^cr(n) with canalizing 
input values a = (ai, . . . , 0,^) and canalized values /3 = (61, . . . , bn)- 
Given a vector a = (oi, 02, . . . , On)-, we define 

OL . . . , J - - - 1 Oiii, , . . . , On ) • 

Then, from the above definition, we immediately have the following result. 

Proposition 3.4. The function f is {a : a : (3} NCF if and only if f is 
{a : a" : NCF. 

Example 3.5. The function fi{xi, X2,X3) = Xi{x2 l)x3 (B 1 is {(1,2,3): 
(0, 1, 0) : (1, 1, 1)} NCF. Actually, one can check that this function is nested 
canalizing in any variable order. Its truth table is given in TableUl 

Example 3.6. Let f2{xi,X2,X3) = {xi(Bl){x2{xsQl)(Bl)®l. This function 
is {(1, 2, 3) : (1, 0, 1) : (1,0, 0)} NCF. It is also {(1,3, 2) : (1, 1, 1) : (1,0, 1)} 
NCF. One can check that this function can be nested canalizing in only two 
variable orders, namely (xi, X2, x^) and (xi, X3, X2). See its truth table in 
Tablel^ 
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Xi 


X2 


X3 


fl 











1 








1 


1 





1 





1 





1 


1 


1 


1 








1 


1 





1 





1 


1 





1 


1 


1 


1 


1 



Xi X2 X3 








1 
1 
1 
1 



Table 1. Truth table for Example 13.51 (left) and Example 13.61 (right) 



Prom the above definitions, we know that, if a function is NCF, all the 
variables appearing in it must be essential. However, a constant function b 
can be < i : a : 5 > NCF for any i and a. 



4. A Detailed Categorization of Nested Canalizing Functions 

As we will see, in a nested canalizing function, some variables are more 
dominant than others. We will classify all the variables of an NCF into 
different levels according to the extent of their dominance. 

Definition 4.1. [18j A function M{xi, . . . ,x„) is an extended monomial of 
essential variables xi, . . . ,Xn if 

M{xi, . . . , Xn) = (xi e ai) • • • {xn © a„), 

where ai S F2 for i G {1, . . . , n}. 

We will rewrite Theorem 3.1 in |17j with more information in the main 
theorem of this section. Basically, we will obtain a unique (the old one is 
not) algebraic normal form (polynomial form). Because of the uniqueness, 
the enumeration of the number of nested canalizing functions, computation 
of their Hamming weight, as well as activity and average sensitivity can be 
done. Besides, in the old form, the variable order of a nested canalizing 
function is not unique. 

Lemma 4.2. The function /(xi, X2, . . . , x„) is < i : a : b > canalizing if 
and only if f{xi, . . . ,x„) = (xj0a)(5(xi, . . . , Xj_i, Xj+i, . . . ,x„)©6 for some 
polynomial Q. 

Proof. Using the algebraic normal form of /, we rewrite it as 



/ — Xigi{xi, . . . , Xj_i, Xj-i-i, . . . , Xn) ® go{xi, ■ ■ ■ , Xj_i, Xj-i-i, . . . , Xn), 
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where gi and go are the quotient and remainder of / when divided by Xj. 
Hence, 

f{X) = f{xi, . . . = (xj e a)gi{xi, . . . . . . 

®agi{xi, . . . ,Xi-i,Xi+i, ... ,Xn) 

®go{xi, . . . ,Xi-i,Xi+i, . . . ,Xn). 

Let gi{xi, . . . ,Xi-i,Xi+i, . . . ,x„) = Q{xi, . . . ,Xj_.i,Xi+i . . . and 

r(xi , . . . , , , . . . , Xn) — Cigi (xx , . . . , Xj_x , Xj^i ; • • • ) 

® go{xi, . . . ,Xi-i,Xi+i, . . . ,Xn)- 

Then 

f{X) = /(xi, ...,Xn) ={xi e a)Q{xi, Xi_i, Xj+i . . . , x„) 

ffi r(xi, . . . , Xj_i, Xj-|_i, . . . , Xn)- 

Since f{X) is < i : a : b > canahzing, we get that 

/(xi, . . . ,Xi_i,a,Xj+i, . . . ,Xn) = b 

for any xi, . . . ,Xj_i,Xi+i, . . . ,Xn, i.e., r(xi, . . . ,Xi_i,Xj+i, . . . ,x„) = 6 for 
any xi, . . . ,Xj_i,Xi+i, . . . ,x„. So r(xi, . . . ,Xi_i,Xi+i, . . . ,x„) must be the 
constant b. This shows necessity, and sufficiency is obvious. □ 

Remark 4.3. We have the following observations. 

(1) When we contrast this lemma to the first part of Theorem 3.1 in |17j . 
it is important to note that Xj is not essential in Q. 

(2) In [30J, there is a general version of this result over any finite field. 

(3) In the above lemma, if f is constant, then Q = 0. 

From Definition 13.31 we have the following result. 

Proposition 4.4. Let /(xi, . . . ,x„) be {a : a : fi} NCF, i.e., f is NCF in 
the variable order Xo.(i) , . . . , 2;o-(n) ^^^^ canalizing input values q = (ai , . . . , On) 
and canalized output values (3 = . . . , 6„). For 1 < k < n — 1, let 

Xa{l) = OT, . . . , Xo-(fc) = Ok- 

Then the function /(xi, . . . , ai , . . . , , . . . , x„) is jcr* : a* : f3*} NCF on 
the remaining variables, where a* = x^(^k+i)j • • • ; x^(n)7 ct* = (ofe+i, . . . , On) 
anc? (3* = (fefe+i, . . ■,bn). 

Definition 4.5. ie^ /(xi,...,Xn) be NCF. We call a variable Xi a most 
dominant variable of f if there is a variable order a = (xj, . . .) such that f 
is NCF in this variable order. 

In Example 13. 5^ all three variables are most dominant. In Example 13. 6^ 
only xi is a most dominant variable. We have: 

Theorem 4.1. Given an NCF /(xi, . . . ,x„), all variables are most domi- 
nant if and only if f = M(xi, . . . , x„)©6, where M is an extended monomial, 
i.e., M = (xi © ai)(x2 ® 02) • • • (x„ © a„). 
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Proof. If xi is most dominant, from Lemma I4.2| we know there exist ai 
and b such that f{xi,X2, • • • , Xn) = {xi © ai)Q{x2, . . . , x„) ® 6, i.e., {xi ® 
ffi Now, if rE2 is also most dominant, then there exist 02 and b' 
such that /(xi, 02, xs, . . . , x„) = 5' for any xi, X3, . . . , x^. Specifically, if 
xi = ai, we get /(ai, 02, X3, . . . , x„) = b = b'. Hence, we also get (x2 ® 
fl2)|(/©^) = {xi (B ai)Q{x2, . . . ,x„). Since xi ©ai and X2 ©02 are coprime, 
we get {x2 ® a2)\Q{x2, . . . ,Xn), hence, /(xi, X2, . . . , x„) = (xi © ai)(x2 © 
a2)Q'(x3, . . . , Xn) © ^. Using the induction principle, the necessity is proved. 
The sufficiency if evident. □ 

We are now ready to prove the main result of this section. Basically, 
we will obtain a new polynomial form of a nested canalizing function by 
induction. In this form, all the variables will be classified into different 
layers, with the variables in the outer layers more dominant than those in the 
inner layers. Variables in the same layer have the same level of dominance. 
Each layer is an extended monomial of the corresponding variables. 

Theorem 4.2. Given n > 2, the function /(xi, . . . , x„) is nested canalizing 
if and only if it can be uniquely written as 

/(xi, . . . , x„) = Mi(M2(- • • (M,_i(M, © 1) © 1) • • • ) e 1) © ^ (4.1) 

where each Mi is an extended monomial. For i ^ j, their corresponding sets 
of variables are disjoint. More precisely, Mi = Y['j=i{^ij © Q^ij ); ^ = !> • • • > 

> 1 for i = 1, . . . ,r - 1, kr > 2, ki + ■ ■ ■ + kr = n, ai^ G F2, {ij \j = 
1,. . . ,ki,i = 1,. . . ,r} = {1,. . . ,n}. 

Proof. We use induction on n. When n = 2, there are 16 Boolean functions, 
8 of which are NCFs, namely 

(xi © ai)(x2 © 02) © c = Ml © 1 © 6, 

where 6 = 1 © c and Mi = (xi © ai)(x2 ffi 02). 

If (xi © ai)(x2 © 02) © c = (xi © ai'){x2 © 02') © c', then, by equating 
coefficients, we immediately obtain ai = ai' , 02 = 02' and c = c'. So 
uniqueness holds. We have proved that Equation 14 . 1 1 holds for n = 2, where 
r = 1. 

Assume now that Equation 14.11 is true for any nested canalizing function 
which has at most n — 1 essential variables. Consider a nested canaliz- 
ing function /(xi, . . . , x„). Suppose Xo-(i) , . . . , Xo-(fc^) are all most dominant 
canalizing variables of /, for 1 < ki < n. 

Case 1: ki = n. Then, by Theorem 14. 11 the conclusion is true with r = 1. 

Case 2: ki < n. Then, with the same arguments as in Theorem 14. 1^ we 
get that / = Mi5r © b, where Mi = (x^(i) © a^(i)) • • • {x^(ki) © ao-(fci))- Let 
Xa{i) = . . . , x^(fcj) = a^(fci) in /, then f = g@b (hence, g) wih also be 
nested canalizing in the remaining variables, by Proposition 14.41 

Since g has n — ki < n — 1 variables, by the induction assumption, we 
get that g = M2(M3(- • • {Mr-i{Mr © 1) © 1) • • • ) © 1) © ^'i- It follows that 
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bi must be 1. Otherwise, all the variables in M2 will also be most dominant 
variables of /. This completes the proof. □ 

Because each nested canalizing function can be uniquely written in the 
form im and the number r is uniquely determined by /, we can make the 
following definition. 

Definition 4.6. For a nested canalizing function f, written in the form 
\4-l[ the number r will be called its layer number. Essential variables of Mi 
will be called most dominant variables (canalizing variables), and are part 
of the first layer of f . Essential variables of M2 will be called second most 
dominant variables and are part of the second layer; etc. 

The function in Example 13.51 has layer number 1, and the function in 
Example 13.61 has layer number 2. 

Remark 4.7. We make some remarks on Theorem \4-.S\ 

(1) It is impossible that kr = 1. Otherwise, Mr © 1 will be a factor of 
Mr^i, which means that the layer number is r — 1. Hence, kr >2. 

(2) If variable xi is in the first layer, and Xi © Oj is a factor of Mi, then 
this nested canalizing function is < i : ai : b > canalizing. We simply 
say that xi is a canalizing variable. 

From the previous examples, we know that a function can be nested canal- 
izing for different variable orders, but only the variables in the same layer 
can be reordered. More precisely, we have the following result. 

Corollary 4.8. If a and t] are two permutations on {l,...,n}, and f is 

both {a : a: (3} NCF and {rj : a' : /3'} NCF, then we have 

{2^(7(1)) • • • ^Xa{ki)} = • • • ;2;^{A:i)}) • • • ; 
{^ct(A;iH hfcr_i+l)' • • • 1 ^o-(ra)} = {2^r)(feiH hfe^-i+l)' • • • ' ^»7{n)}- 

We also have /3 = /3' = (61, . . . , and 

bi = ■ ■ ■ = bk^, . . . ,6a:iH hA:,.-i+l = ' ' ' = &n- 

Furthermore, 

bl / fefci+l / fefci+fca+l / • • • / bki+-+kr+i+l- 
Proof. These results all follow from the expression for / in Theorem 14.21 □ 

From this corollary we can determine the layer number of any nested 
canalizing function by its canalized value. For example, if / is nested canal- 
izing with canalized value (1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 1) (n = 11), then its layer 
number is 5. 

Let NCF(n, r) denote the set of all nested canalizing functions in n vari- 
ables with layer number r, and let NCF(n) denote the set of all nested 
canalizing functions in n variables. 
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Corollary 4.9. For n >2, 

|MCF(n,r)| = 2"+i E "fe J' 



fciH \-kr=n 

ki>l,i=l,...,r-l,kr>2 



and 

n-1 



incf(„)i=2»«j: ^ u,„,,J 

r=l ki+---+kr=n ^ ' »- J-/ 

A;i>l,j=l,...,r-l,fcr>2 

where the multinomial coefficient (^^ ^) is eguaZ to ^^"'^ , . 

Proof. It follows from Equation 14.11 that for each choice of ki, . . . ,kr, with 
A;i+- • •+A;r = n, fe^ > 1, i = 1, . . . , r-1 and kr > 2, there are 2^^ ^n-k^-^-k,.^^ 
ways to form Mj, j = 1, . . . ,r. 

Note that we have two choices for b. Hence, 



|NCF(n,r)| = 2 ^ 



■+kr 



fciH \-kr=n 

ki>l,i=l,...,r-l,kr>2 

n\ fn — ki\ /n — ki — ■ ■ ■ — kr-i 
hj \ h J \ kr 

= 2"+i Yl 

fclH \-kr=n 

k^>l,i=l,...,r-l,kr>2 

n\ {n — ki)\ 

{kiy.{n-kiy.{k2)l{n-ki-k2)l"' 

{n — ki — ■ ■ ■ — kr-i)\ 
kr\{n — ki — ■ ■ ■ — krY- 



ki>l,i=l,...,r-l,kr>2 



^ \ki, . . . , kj.-i 

ki-^ \-kr=n 

k,>l,i=l,...,r-l,kr>2 



Since NCF(n) = UP=i NCF(n,r) and NCF(n, i) fl NCF(n, j) = </> when i / 
j, we get the formula for |NCF(n)|. □ 

As examples, one can check that |NCF(2)| = 8, |NCF(3)| = 64, |NCF(4)| = 
736, |NCF(5)| = 10624, .... These results are consistent with those in [4ti38j. 
By equating our formula to the recursive relation in [U [38] , we have 
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Corollary 4.10. The solution of the nonlinear recursive sequence 

n-l 

a2 



r=2 

is 



n-l 

n 



r=l fciH \-kr=n 

fc>l,i=l,...,r--l,fcr>2 



/ci , . . . , kr—l 



5. Hamming Weight, Activity, and Average Sensitivity 

A Boolean function is called balanced if it takes the value 1 on exactly half 
the states (and on the other half). In other words, its Hamming weight (the 
number of I's in its truth table) is 2*^"^, where n is the number of variables. 
Hence, there are (2^-1) balanced Boolean functions. It is easy to show that 
a Boolean function with canalizing variables is not balanced, i.e., is biased, 
actually, very biased. For example, the two constant functions are trivially 
canalizing, and they are also the most biased functions. Extended monomial 
functions are the second most biased since for any of them, only one value 
is nonzero. But biased functions may have no canalizing variables. For 
example, /(xi, 2:2, ^3) = xiX2X^ © a:;ia;2 © xix^ © X2X^ is biased but without 
canalizing variables. 

In Boolean functions, some variables have greater influence over the out- 
put of the function than other variables. To formalize this, a concept called 
activity was introduced. Let x = (rci, . . . ,x„), and 

— TT — ^ = /(s^i, • • • 5 © 1) • • • ) 3;„) © /(xi, . . . , Xj, . . . , x„). 

OXi 

The activity of variable Xj is defined as 

1 9/(xi,...,Xn) 

(xi,...,x„)GF5 

Note that the above definition can also be written as follows: 
H = E (/(xi, . . . ,0, . . . ,x„)©/(xi, . . . , 1, . . . ,Xn)) 



(xi,...,a;i_i,a;i+i,...,a;„)GF2 



(5.2) 

The activity of any variable in a constant function is 0. For an affine function 
/(xi, . . . , Xn) = xi © ■ • ■ © Xn © 6, A{ = 1 for any i. It is clear, for any / and 
i, that we have < A{ < 1. 

Another important quantity is the sensitivity of a Boolean function, which 
measures how sensitive the output of the function is if the input changes 
(This was introduced in p^). The sensitivity s-'^(xi, . . . , x„) of / on the in- 
put (xi, . . . , Xn) is defined as the number of Hamming neighbors of (xi, . . . , x^) 
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(that is, all states that have Hamming distance 1) on which the function 
value is different from f{xi, . . . ,Xn)- That is, 

S^{xi,...,Xn) = \{i\f{xi,...,0,...,Xn) + 

i 

f{xi,.. .,Xn),i = 1, • • ■,ri}\. 

Obviously, s/(xi, ...,xn) = Er=i ^%f^- 

The average sensitivity of a function / is defined as 

= E[s^{xi, . . . ,x„)] = 
1 " 

— Yl sf{xi,...,Xn)=Y,^(- 

(xi,..,x„)GF^ i=l 

It is clear that < s-'^ < n. The concept of average sensitivity of a Boolean 
function is one of the most studied concepts in the analysis of Boolean 
functions, and has received a lot attention recently [H [SI El [71 [lOl [HI [Ml 
[26l[28l[37l[39l[40l[4ll[42l[44]. Bernasconi has shown [5] that a random 
Boolean function has average sensitivity ^. This means the average value 
of the average sensitivities of all Boolean functions in n variables is ^. In 
|40j , Shmulevich and Kauffman calculated the activity of all the variables of 
a Boolean function with exactly one canalizing variable and unbiased input 
for the other variables. Adding all the activities, the average sensitivity of 
a Boolean function was also obtained. 

First, the following observation will be useful. We have the equality 

1, if (xi, . . . jXfc) = (oT, . . . ,a^) 
0, otherwise. 

That is, only one value is equal to 1 and all the other 2^ — 1 values are 0. 



(xi e ai) • • • (xfc e ofc) 



Theorem 5.1. For n>2, let /i = Mi, 

fr = Mi(M2(- • • {Mr-l{Mr 1) 1) • • • ) 1), r > 2 
where Mi is same as in Theorem \4-^ Then the Hamming weight of fr is 



W{fr) = ^(-l)-'-i2"-n=i ''^ (5.3) 
The Hamming weight of fr (Bl is 



r 



W{fr e 1) = ^(-l)^^2"-S-i'=% (5.4) 
where Y^=i ^-^ interpreted as 0. 

Proof. First, consider the Hamming weight of fr- When r = 1, we know the 
result is true by the above observation. When r > 1, we consider two cases: 
Case ^: r = 2i + 1 is odd. Then all the states on which / evaluates to 1 
will be divided into the following disjoint groups: 
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• Group j for j = 1, . . . , t : Ml = 1, M2 = 1, . . . , M2j^i = 1, M2j = 0; 

• Group t + 1 : Ml = 1, M2 = 1, . . . , M2t = 1, M2t+i = Mr = 1. 
In Group j, the number of states is 

^2^2j _ \'j2"'^^'^ '^^j — 2"-~'i^l ^2j-l _ r^—kx kij 

In Group t + the number of states is 2^^^^ = 1. 

Adding all of them together, we get Equation 15.31 

Case B: r = 2t is even. The proof in this case is similar, and we omit it. 
Because |{(xi, . . .,Xn)\f{xi, . . . ,x„) = 0}| + |{(xi, . . . ,a;„)|/(xi, . . . ,x„) = 
1}| = 2", we know that the Hamming weight of © 1 is equal to 

r 

W{fr © 1) = 2" - W{fr) = 2" - ^(-l)i-l2"-^-i'=' 

i=i 

r 

= ^(-l)J2"-^-i'=% 

j=0 

where Yl'i=i should be interpreted as 0. □ 

In the following, we will calculate the activities of the variables of any 
nested canalizing function. For this, we will use the formula for the Hamming 
weight of a nested canalizing function since the function in the summation 
will be reduced to a nested canalizing function, for which the first layer is 
a product of a few layers of the original nested canalizing function. Let 
/ be nested canalizing, written in the form of Theorem 14.21 Without loss 
of generality, we can assume that Mi = (xi © ai)(x2 © 02) • • • (xk^ © a^^). 
Let mi = {xi © ai) ■ ■ ■ (xj_i © aj_i)(j;j+i © Oj+i) • • • {xk-^ © a^J, so that 
Ml = {xi © ai)mi. 

If r = 1, i.e., ki = n, then 

{xi,...,Xi-i,Xi+i,...,x„)e¥2~^ 

i i 
(/(Xi, . . . , 0, . . . , X„) © /(Xi, Xn)) 

1 



> mi 



2n-i 

(x-l,...,Xi-l,Xi + i,...,Xn)&2^^ 

-}—W(mi) = 

2?!— 1 ^ ' 2"~ 

If 1 < r < n — 1, then consider the activity of xi in the first layer, i.e., 
1 < i < ki. We have 

= ^— V 

i 2"-l ^ 

{x-x,...,Xi-l,Xi+i,...,Xn)&2~^ 

i i 
(/(Xi, . . . , 0, . . . , Xn) © /(Xi, . . . , 1, . . . , Xn)) 
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1 



2"— 1 ^ ^ 

mi(M2(- • • {Mr-i{Mr e 1) e 1) • • • ) e 1) 



{xi,...,Xi-l,Xi+-i,...,Xn)&2 



-L^wimi{M2{- ■ ■ (M,_i(M, e 1) e 1) • • • ) e i)). 





-1) 


^f 


ki > 1 






if 


ki = 1 


f 2^E,^=i(-ir'2«-n^i'=% 


^/ 


ki 


> 1 


I 2^E,^=o(-ip2"-^-»'-s 


^/ 


ki 


= 1. 


/ 2^E,'=i(-ir'2«^n^i'=s 


i/ 


ki 


> 1 


1 2^E,^=^(-iy2"-^-s 


^/ 


ki 


= 1. 


/ 2^E-=i(-ir'2«-n^i'=% 


if 


ki 


> 1 


1 2^E,'=i(-ir'2«-n^i'=s 


if 


ki 


= 1. 



= 2^Ej=i(-l)^'~^2"-S,=ife> by TheoremEU Note, in the above, ki = 
1 means mi = 1, so we used Equation 15. 41 with layer number r — 1, and the 
first layer is M2 for n — 1 variable functions. 

Now let us consider the variables in the / — th layer, i.e., Xi is an essential 
variable of M;, 2 < / < r — 1. We have Mi = (xj + ai)mi and 

{xi,...,Xi-l,Xi+i,...,Xri)e¥2~^ 

i i 
(/(Xi, . . . , 0, . . . , X„) e /(Xi, . . . , 1, . . . , Xn)) 

= ^ y 

2"— 1 / > 

{xi,...,Xi-l,Xi+i,...,Xn)eF2~^ 

Ml • • • Mi_imi{Mi+i{- • • (M^ e 1) • • • ) e 1). 

, r-l+l 

= ^ (_X)i-i2""^"((^i+'"+'''"^)+*='+i+'"+''j+'-i)) = 

2n— 1 / ^ ^ ' 

i=i 

by Equation 15.31 in Theorem 15. li Note that Mi ■ ■ ■ Mi-iuii is the first layer, 
M;_|_i is the second layer, etc. 

Let Xi be the variable in the last layer M^, then we have 

^ ^ M1M2 ■ ■ ■ Mr-irur ^ 



2n-l ^.^i^.^z 2"-l' 

(xi,...,Xi_i,Xi+i,...,a;„)GF2~"'" 

Variables in the same layer have the same activities, so we use Aj to stand 
for the activity number of each variable in the Ith layer Mi, 1 < I < r. We 
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find that the formula of Aj for 2 < / < r — 1 is also true when / = r or r = 1. 
The next theorem summarizes these. 

Theorem 5.2. Let f be a nested canalizing function, written as in Theorem 
\4-^ Then the activity of each variable in the Ith layer , 1 < / < r, is 

r-l+l 

M = i^iY. (5.5) 

The average sensitivity of f is 

sf = ±kA{ = ^E^' EVlr^2"-E£^^'=^ (5.6) 
1=1 1=1 j=i 

Next, we analyze the formulas in Theorem 15.21 

Corollary 5.1. If n > 3, then A( > > ■ ■ ■ > aI , and < s^ < 

9 1 

Proof. 

r-l+l 



2 



^_(2"~^^ ^' — 2"~^i ^'+1 + • • • (—1)''^') 



2^^ 

Since the sum is an alternating decreasing sequence and ki^i > 1, we have 
^ /-pw-fci ki-l\ ^ ^ /2"~fei '^i — 2"~^i ^i+i) 

on— 1 ^ ' — 9?i— 1 ^ ' 



2*^ 

I 9?i— 1 ^ ' 



Hence, 



< J_(2"-'=i— < a{ 

— 9n— 1 ^ / ^ i 



n-fci fci+i^ 

2*^ 



We have 



k^A{ = ^^(2"-'=i - 2"-^i-'=2 + . . . (-1)^-1); 



^^^/ _ ki _ 2"-^! ki~ki+i _ . . . (^_iY~iy 



2^^ 



2n-l- 

if J^l 



Hence, = YlUi ^iM ^ 2^ + 2^ + • • • + 2^ = so we know 

that the nested canalizing functions with layer number 1 have minimal 
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average sensitivity. On the other hand, = ^^1=1 ^ 2" ^'^ + 
+ . . . + ^2''-k,--h + . . . + = u{ku ■■■ , kr), where 
ki + ■ ■ ■ + kr = n, ki > 1, i = 1, . . . , r — 1 and k^ > 2. We wih find the 
maximal value of U{ki, . . . , kr) in the following. 

First, we claim /c^ = 2 if U{ki, . . . , kr) reaches its maximal value. Because, 
if kr is increased by 1, and the last term makes more contributions 
to U{ki, . . . ,kr), then there exists / such that ki will be decreased by 1 
{ki + ■ ■ ■ + kr = n), hence 

2"~^'i h 

will be decreased more than 2!r~r- Now, it is obvious that ^^2^~''^ attains 
its maximal value only when ki = 1 01 2, but ki = 1 will be the choice since 
it also makes all the other terms greater. Likewise, ^^li 2"~^i~^^ attains 
its maximal value when ki = k2 = 1 or ki = 1 and k2 = 2; again, /c2 = 1 
is the best choice to make all the other terms greater. In general, if ki = 
■ ■ ■ = ki_i = 1, then ^^1^ 2""^^ ^' attains its maximal value when ki = 1, 
where 1 < / < r — 1. In summary, we have shown that U (/ci, . . . , kr) reaches 
its maximal value when r = n — 1, ki = ■ ■ ■ = kn-2 = 1, kn-i = 2, and 

Max[U{ki,. ..,kr) = U{1, . . . , 1, 2)] = 
1 /2"-i + 2"-2 + . . . + 2^ + 2) = 2 - 

2n-l ^ > > > J 2"~^ 

□ 

Remark 5.2. The minimal value of average sensitivity approaches and the 
maximal value of U{ki, . . . , kr) approaches 2 as n — > 00. Hence, < < 2 
for any NCF with an arbitrary number of variables. 

In the following, we evaluate Equation 15.61 for some parameters ki, . . . ,kr. 

Lemma 5.3. (1) If r = n — l,ki = ■■■ = kn-2 = Ij^n-i = 2, then 

/ _ 4 _ 3+(-l)" 
* ~ 3 3x2" • 

(2) Given n > 4, r = n — 2,ki = ... = kn-3 = l,kn-2 = 3, then 

/ _ 4 _ 9+5(-l)"-i 
^ 3 3x2" 

(3) Ifn is even and n > 6, r = ^, ki = 1, k2 = • • • = ^^-1 = 2, /c^ = 3, 
then s^ = ^ — j^2^ . Hence, these three cardinalities are equal if n is 
even. 

Proof. When r = n — 1, then ki = ■ ■ ■ = kn-2 = ^n-i = 2 by Equation 
15. 6[ We have 

r n—l n—l 

1=1 1=1 j=i 

= ^E^'( E (-ir^2"-^-'+i + (-ir-'-i) 
1=1 j=i 
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1=1 

= J_(y (i2"-'+i + -(-1)"-') + 2) = - - ^lizDl. 

2n-i^Z^V3 ^\ J J ; 3 3^2" 

1=1 

The other two formulas are also routine applications of Equation 15.61 □ 

Based on our numerical calculations, Lemma |5.3| and the proof of Corol- 
lary [5TT1 we can make the following conjecture. 

r f A 3+( 1)^ 

Conjecture 5.4. The maximal value of is = 3x2" • 

he reached if the nested canalizing function has the maximal layer number 
n — 1, i.e., if r = n — \, ki = ■ ■ ■ = kn-2 = 1; ^n-i = 2. When n is 
even, this maximal value is also reached by a nested canalizing function with 
parameters n > A, r = n — 2, ki = ■ ■ ■ = kn-3 = 1, kn-2 = 3 or n > 6, 
r = ^ , ki = 1, k2 = ■ ■ ■ = k^_i = 2 and kii = 3. 

Remark 5.5. When n = 6, the nested canalizing function with ki = 1, 
k2 = 2, k^ = 1, and k4 = 2 also has the maximal average sensitivity 
But this can not be generalized. If the above conjecture is true, then we have 
< s-^ < I for any nested canalizing function with an arbitrary number of 
variables. In other words, both and | are uniform tight bounds for any 
nested canalizing function. 

6. Simulations of Network Dynamics 

The sensitivity of functions [40j has been shown to be a good indicator 
of the stability of dynamical networks constructed using random Boolean 
functions. We have generated random networks controlled by nested canal- 
izing functions with fixed layer number m, where m = 1, . . . ,n — 1. For 
our simulations we followed a similar approach as in [30]. Starting with a 
Boolean network F, we sample pairs of random states x and y for F. Let 
H(t) be the Hamming distance of x and y, i.e. II(t) is the number of bits 
in which x and y differ, and let H{t + 1) be the Hamming distance of F(x) 
and F{y), i.e. the Hamming distance of the successor states of x and y. A 
Derrida curve [13] is a plot of H{t + 1) against II{t) for all possible Ham- 
ming distances. Figure [1] shows Derrida plots for different layer numbers. 
Each Derrida curve was generated from 4096 random networks by taking 
the average Hamming distances. 

As can be observed from Figure [H networks made up of nested canaliz- 
ing functions of fixed layer number equal to 1 are significantly more stable 
than networks constructed from nested canalizing functions with higher layer 
numbers. It should be noted that the class of nested canalizing functions 
with layer number equal to 1 is the same family as the AND-NOT networks 
studied in [46i, equal to the class of extended monomials. Similarly, net- 
works made up of nested canalizing functions with fixed layer number equal 
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Derrida plot for NCFs with n=100 and k=12 




Initial Hamming distance 

Figure 1. Derrida plots for different layer numbers. Here, 
n is the number of nodes and k is the connectivity. For 
each layer number, a Derrida plot was generated from 4096 
random networks. The x-axis represents Hamming distance 
of pairs of states, and the y-axis represents Hamming distance 
of their images. 

to 2 are more stable than networks constructed from nested canalizing func- 
tions with higher layer numbers, and so on. Figured] shows that, as the layer 
number increases, networks becomes less stable. This matches our results 
in Section [5l 

7. Conclusion 

Nested canalizing Boolean functions were inspired by structural and dy- 
namic features of biological networks. In this study we took a careful look 
at the computational properties of Boolean functions and of Boolean net- 
works constructed from them. The main tool for our analysis is a particular 
polynomial normal form of the Boolean functions in question. In particular, 
we introduced a new invariant for nested canalizing functions, their layer 
number. Using it, we obtain an explicit formula of the number of nested 
canalizing functions, which improves on the known recursive formula. Based 
on the polynomial form, we also obtain a formula for the Hamming weight 
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of a nested canalizing function. The activity number of each variable of 
a nested canalizing function is also provided with an explicit formula. An 
important result we obtain is that the average sensitivity of any nested 
canalizing function is less than 2, which provides a theoretical argument 
why nested canalizing functions exhibit robust dynamic properties. This 
leads us to conjecture that the tight upper bound for the average sensitivity 
of any nested canalizing function is |. 

It should be noted that all the variables in the first layer of an nested 
canalizing function are canalizing variables (the most dominant), all the 
variables in the second layer arc the second most dominant, etc. The fact 
that networks with low layer number are more stable than those with high 
layer number could be a consequence of this observation, since lower layer 
number means more dominant variables. The most extreme examples are: 
when the layer number is equal to 1, we have proved that such a function 
has minimal average sensitivity, therefore we also conjecture that the func- 
tion with layer number number equal to tt, — 1 (the maximal layer number) 
has maximal average sensitivity. Hence, it should be reasonable that the 
corresponding networks show similar behavior. 

Acknowledgments 

The authors thank the referees for insightful comments that have im- 
proved the manuscript. 

References 

[1] Kazuyuki Amano, "Tight bounds on the average sensitivity of k-CNF, "Theory of 
Computing, Vol 7 (2011), pp. 45-48. 

[2] E. Balleza, E. R. Alvarez-Buylla, A. Chaos, S. Kauffman, I. Shmulevich,and M. Al- 
dana, "Critical dynamics in genetic regulatory networks: Examples from four king- 
doms, "PLoS ONE, 3 (2008), p. e2456 

[3] C. Barrett, C. Herring, J. Reed, and B. Palsson, "The global transcriptional regu- 
latory network for metabolism in Escherichia coli exhibits few dominant functional 
states, "Proc Natl Acad Sci USA, 102 (2005), pp. 1910319108. 

[4] E. A. Bender, J. T. Butler, "Asymptotic approximations for the number of fanout-tree 
functions, "IEEE Trans. Comput. 27 (12) (1978) 1180-1183. 

[5] A. Bernasconi, "Mathematical techniques for the analysis of Boolean functions, 
"Ph.D. thesis, Dipartmento di Informatica, Universita di Pisa (March, 1998). 

[6] A. Bernasconi, "Sensitivity vs. block sensitivity (an average-case study) "Information 
processing letters 59 (1996) 151-157. 

[7] Ravi B. Boppana, "The average sensitivity of bounded-depth circuits "Information 
processing letters 63 (1997) 257-261. 

[8] J. T. Butler, T. Sasao, and M. Matsuura, "Average path length of binary decision 
diagrams, "IEEE Transactions on Computers, 54 (2005), pp. 10411053. 

[9] David Canright, Sugata Gangopadhyay, Subhamoy Maitra, Pantelimon Stonica, 
"Laced Boolean functions and subset sum problems in finite fields, "Discrete applied 
mathematics, 159 (2011), pp. 1059-1069. 
[10] Shijian Chen and Yiguang Hong, "Control of random Boolean networks via average 

sensitivity of Boolean functions, "Chin. Phys. B Vol. 20, No 3 (2011) 036401. 
[11] Demetres Christofides, "Influences of Monotone Boolean EXinctions, "Preprint 2009. 



20 Y. LI, J. O. ADEYEYE, D. MURRUGARRA, B. AGUILAR, R. LAUBENBACHER 

[12] S. A. Cook, C. Dwork, R. Reischuk, "Upper and lower time bounds for parallel 
random access machines without simultaneous writes, "SIAM J. Comput, 15 (1986), 
pp. 87-89. 

[13] B. Derrida and Y. Pomeau, "Random networks of automata: a simple annealed 

approximation", Europhys, Lett., 1:45-49, 1986. 
[14] S. E. Harris, B. K. Sawhill, A. Wuensche, and S. Kauffman, "A model of transcrip- 
tional regulatory networks based on biases in the observed regulation rules, " Complex, 
7 (2002), pp. 2340. 

[15] M. Herrgard, B. Lee, V. Portnoy, and B. Palsson, "Integrated analysis of regula- 
tory and metabolic networks reveals novel regulatory mechanisms in saccharomyces 
ccrcvisiac, "Genome Res, 16 (2006), pp. 627635. 

[16] A. Jarrah, R. Laubcnbaclicr, and A. Vcliz-Cuba, "A polynomial framework for mod- 
eling and anaylzing logical models. "In Preparation, 2008. 

[17] A. Jarrah, B. Ropasa and R. Laubenbacher, "Nested canalyzing, Unate Cascade, and 
Polynomial Functions", Physica D 233 (2007), pp. 167-174. 

[18] Winfried Just, "The steady state system problem is NP-hard even for monotone 
quadratic Boolean dynamical systems "Preprint ,2006 

[19] Winfried Just, Ilya Shmulevich, John Konvalina, "The number and probability of 
canalyzing functions", Physica D 197 (2004), pp. 211-221. 

[20] S. A. Kauffman, "The Origins of Order: Self-Organization and Selection in Evolu- 
tion", Oxford University Press, New York, Oxford (1993). 

[21] S. A. Kauffman, C. Peterson, B. samuelesson, C. Troein, "Random Boolean Network 
Models and the Yeast Transcription Network", Proc. Natl. Acad. Sci 100 (25) (2003), 
pp. 14796-14799. 

[22] S. A. Kauffman, C. Peterson, B. Samuelsson, and C. Troein, "Genetic networks with 
canalyzing Boolean rules are always stable, ", PNAS, 101 (2004), pp. 1710217107. 

[23] S. A. Kauffman, "The large-scale structure and dynamics of gene control circuits: an 
ensemble approach, ", J. Theor. Biol. 44 (1974) 167. 

[24] N.Keller and H. Pilpel, "Linear transformations on monotone functions on the discrete 
cube, "Discrete Math. 309 (2009), 4210-4214. 

[25] R. Laubenbacher and B. Pareigis, "Equivenlence relations on finite dynamical sys- 
tems, ", Advances in applied mathematics 26 (2001), pp. 237-251. 

[26] W. Liu, H. Lohdcdmaki, Edward R. Dougherty and I. Shmulevich, "Inference of 
Boolean Networks Using Sensitivity Rcgularization, ", EURASIP Journal of Bioin- 
formatics and System Biology Volume 2008, Article ID 780541, 12 pages. 

[27] Layne, L., Dimitrova, E.S., Macauley, M. (2012). Biologically relevant properties of 
nested canalyzing functions. Bulletin of Mathematical Biology, 74(2), pp. 422-433. 

[28] Jiyou Li, "On the average sensitivity of the weighted sum function, 
"arXiv:1108.3198v2 [cs.IT] 18 Aug 2011. 

[29] Yuan Li, "Results on Rotation Symmetric Polynomials Over GF{p)" , Information 
Sciences 178 (2008), pp. 280-286. 

[30] Yuan Li, David Murrugarra, John O Adeyeye and Reinhard Lauben- 
bacher "Multi-State Canalyzing E\inctions over Finite Fields ", 
http://arxiv.org/pdf/1110.6481vl.pdf 

[31] R. Lidl and H. Niederreiter, "Finite Fields", Cambridge University Press, New York 
(1977). 

[32] A. A. Moreira and L. A. Amaral, "Canalyzing Kauffman networks: Nonergodicity 
and its effect on their critical behavior, ", Phys. Rev. lett. 94 (21) (2005), 218702. 

[33] David Murrugarra and Reinhard Laubenbacher, "Regulatory patterns in molecular 
interaction networks, "(2011), Journal of Theoretical Biology, 288, 66-72. 

[34] David Murrugarra and Reinhard Laubenbacher, "The number of multistatc nested 
canalyzing functions, "(2012), Physica D: Nonlinear Phenomena, 241, 921-938. 



BOOLEAN NESTED CANALIZING FUNCTIONS 



21 



[35] S. Nikolajewaa, M. Friedela, and T. Wilhelm, "Boolean networks with biologically 
relevant rules show ordered behaviorstar, open, "Biosystems, 90 (2007), pp. 4047. 

[36] N. Nisan, "CREW PRAMs and decision tree, "SIAM J. Comput, 20 (6) (1991), PP. 
999-1070. 

[37] Xiaoning Qian and Edward R. Dougherty, "A comparative study on sensitivitys of 

Boolean networks, "978-1-61284-792-4/10 2011 IEEE. 
[38] T. Sasao, K. Kinoshita, "On the number of fanout-tree functions and unate cascade 

functions, "IEEE Trans. Comput. 28 (1) (1979) 66-72. 
[39] Ilya Shmulevich, "Average sensitivity of typical monotone Boolean functions, 

"arXiv;math/0507030vl [math. Co] 1 July 2005. 
[40] Ilya Shmulevich and Stuart A. Kauffman, "Activities and sensitivities in Boolean 

network models, "Physical Review Letters, Vol 93, Number 4 (2004), 048701. 
[41] Igor. E. Shpaxlinski, "Bounds on the Fourier coefficients of the weighted sum function, 

"Information Processing Letters. 103 (2007), 83-87. 
[42] Stoffcn Schobcr and Martin Bosscrt, "Analysis of random Boolean networks using 

the average sensitivity, "arXiv: 0704.0197vl [nlin.CG] 2 Apr 2007. 
[43] C. H. Waddington, "Canalisation of development and the inheritance of acquired 

characters, "Nature, 150 (1942), pp. 563564. 
[44] Madars Virza, "Sensitivity versus block sensitivity of Boolean functions, 

"arXiv;1008.0521v2 [cs.CC]8 Dec 2010. 
[45] Shengyu Zhang, "Note on the average sensitivity of monotone Boolean functions, 

"Preprint 2011. 

[46] A. Veliz-Cuba, K. Buschur, R. Hamershock, A. Kniss, E. Wolff, and R. Laubenbacher, 
"AND-NOT decomposition of boolean network models". Preprint 2012. 

^Department of Mathematics, Winston-Salem State University, NC 
271 10, usa, email: liyu@wssu.edu, ^department of mathematics, wlnston- 
Salem State University, NC 27110, USA, email: adeyeyej@wssu.edu, ^ School 
OF Mathematics, Georgia Institute of Technology, Atlanta, GA 30332- 
0160 USA, email: davidmur@math.gatech.edu, ''Department of Com- 
puter Science, Virginia Tech Blacksburg, VA 24061-0123, USA, email: 

BAGUILAR@VT.EDU, ^ VIRGINIA BlOINFORMATICS INSTITUTE, VIRGINIA TECH, 

Blacksburg, VA 24061, US A and Department of Mathematics, Virginia 
Tech, Blacksburg, VA 24061-0123, USA, email: reinhard@vbi.vt.edu 



